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We investigate nucleosynthesis inside both a gamma-ray burst accretion disk and 
a wind launched from an inner region of the disk using one-dimensional models of 
the disk and wind and a nuclear reaction network. Far from a central black hole, 
the composition of accreting gas is taken to be that of an O-rich layer of a massive 
star before core collapse. We find that the disk consists of five layers characterized 
by dominant elements: ls O, 28 Si, 54 Fe (and 56 Ni), 4 He, and nucleons, and the 
individual layers shift inward with keeping the overall profiles of compositions 



as the accretion rate decreases. 56 Ni are abundantly ejected through the wind 



from the inner region of the disk with the electron fraction ~ 0.5. In addition to 
iron group, elements heavier than Cu, in particular 63 Cu and 64 Zn, are massively 
produced through the wind. Various neutron-rich nuclei can be also produced in 
the wind from neutron-rich regions of the disk, though the estimated yields have 
large uncertainties. 



1. Introduction 

Observational evidences have been accumulated for a connection between 
gamma-ray bursts (GRBs) and supernovae (SNe): association of SN 1998bw 
and GRB 980425 1 and SN 2003dh in afterglow of GRB 030329. 2 A col- 
lapsar model is one of promising scenarios to explain a huge gamma-ray 
production in GRBs and GRB/SN connections. 3,4 During collapse of mas- 
sive stars, stellar material greater than several solar masses falls back on 
a new-born black hole with extremely high accretion rates (< lM Q s _1 ). 5 
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An accretion disk forms around the hole due to the angular momentum 
of the fallback material. 3 ' 4,6 In the context of the collapsar model, jet-like 
explosion driven by neutrino annihilation and nucleosynthesis in the jet has 
been investigated. Although 56 Ni with high velocity (> 0.1c) can be mas- 
sively produced , it is not sufficient for an observed amount in SN 1998bw. 7 
In addition to the production via the jet, massive synthesis of 56 Ni is also 
suggested in winds launched from the accretion disk. 3,8 Neutron-rich nuclei 
may be produced through r-process inside the wind ejected from an inner, 
neutron-rich region of the disk. 8 In the present paper, we examine nucle- 
osynthesis inside a GRB accretion disk and investigate abundance change 
through the wind launched from the disk. 

2. Disk Model and Input Physics 

We construct a steady, axisymmetric model 9,10 of the disk around a black 
hole of mass M with accretion rates up to 10A/ Q s _1 . The black hole mass is 
fixed to be 3 M@ and the viscosity parameter is set to be a v i s = 0.1. Figure 

1 shows the profiles of density, p, (thick lines) and temperature, T, (thin 
lines) in the accretion disk for M = 1 (solid lines), 0.1 (dotted lines), and 
O.OIMqS -1 (dashed lines). These profiles arc roughly agreement with the 
corresponding profiles in another disk model 9 except for the region near the 
disk inner edge; our density and temperature drop rapidly near the inner 
edge. This is attributed to our use of the pseudo-Newtonian potential and 
the zero torque condition at the inner edge. 

Chemical composition of accreting gas far from the black hole is set to 
be that of an O-rich layer of a 20 M Q star before core collapse. 11 Once 
temperatures and densities are evaluated inside the disk, using a nuclear 
reaction network, 12,13 which includes 463 nuclei up to Kr, we follow evo- 
lution of the composition in the accreting material during the infall onto 
the hole. We note that at the inner region of the disk, where T > 9 x 10 9 K, 
the composition is in nuclear statistical equilibrium (NSE). 14 

3. Abundance Distribution inside Disks 

Figure 2 shows the abundance profiles of representative nuclei inside the 
disk with M = O.lMgS -1 . Far from the black hole r > 1000r g , where 
r g = 2GM/c 2 is the Schwarzschild radius, accreting gas keeps presupernova 
composition, or that of the O-rich layer, because of low temperatures (T < 

2 x 10 9 K). As the material falls down, the gas becomes rich in iron-group 
elements via explosive O-burning, followed by Si-burning. The processed 
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Figure 1. Density and temperature profiles inside the disk with 0.01, 0.1 and IMq s . 

heavy elements, however, are destroyed to helium, and finally to protons 
and neutrons through photodisintegrations deep inside the disk. Near the 
inner edge of the disk, neutrons are the most abundant by electron captures 
on protons. Considerable amounts of D, T and 6 Li exist due to NSE. 

For higher accretion rates, the density is higher, the electron capture is 
more efficient, and consequently the disk becomes more neutron-rich. In 
fact, the ratios of neutron to proton are 1.33, 3.50 and 10.4 for M = 0.01, 
0.1 and 1 MqS" 1 , respectively, near the inner edge of the disk. Radial 
profiles of the electron fraction are similar to those in the other authors. 8 

The disk consists of five layers characterized by dominant elements: 16 0, 
28 Si, 54 Fe (and 56 Ni), 4 He, and nucleons. Hereafter, these five layers are 
referred to as the O-rich, Si-rich, Fe-rich, He-rich, and np-rich disk layers. 
Figure 3 shows the interfaces of these layers for M = 0.01— IMq s _1 . Tem- 
peratures increase with the increasing accretion rate at a given radius, so 
that the individual layers shift outward, but the overall profiles of compo- 
sition are preserved. This is because nucleosynthesis inside a disk mainly 
depends on the temperature distribution of the disk. Averaged mass frac- 
tions over the individual layers of 40 abundant nuclei are given in Table 
1 for the Fe-rich and He-rich disk layers with M = 0.01 and 0.1 Mq s . 
We find that the averaged abundances of the individual disk layers are not 
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Figure 2. Abundance profiles of representative nuclei inside the disk with O.IMq s 1 . 



significantly changed as the accretion rates. 



4. Chemical Composition of Disk Winds 
4.1. 56 Ni in Winds 

During accretion onto a black hole, some fractions of accreting gas could be 
ejected through winds from an accretion disk. Possible processes driving the 
wind are magnctical centrifugal force 15 and viscosity. 3 Abundances of ejecta 
through the wind change via decay processes as well as charged particle and 
capture processes because of high densities and temperatures of the ejecta. 
The abundances hence depend not only on initial conditions of the ejecta 
(abundances, density, temperature, and so on) but also on hydrodynamics 
of the wind. Detailed dynamics of the wind is, however, still uncertain. 
We therefore adopt a simple hydrodynamical model of the wind, 16 where 
the gasses are assumed to be adiabatic and freely expanding. The ejection 
velocity w c j is set to be 10~ 4 — O.Iuk with the Keplerian velocity vk at 
the ejection radius r j of the disk. The density and temperature of the 
wind are initially taken to be those of the disk at r e j. We note that the 
entropy per baryon are 10-20 in units of the Boltzmann constant and that 
the non-relativistic entropies are always larger than the relativistic ones in 
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Figure 3. Interfaces of various layers . The disk consists of five layers characterized by 
dominant elements: ls O, 28 Si, 54 Fe (and 56 Ni), 4 He, and nucleons. 

the ejecta. The adiabatic index is accordingly taken to be 5/3. 

Using density and temperature evolution calculated with the above wind 
model, we can evaluate change in abundances of the ejecta from an initial 
composition which is the same as in the accretion disk at r e j, with similar 
post-processing calculations to those in accretion disks. We have calculated 
abundance change through wind launched from the inner region (r cj < 
200r s ) of the disk with M = 1CP 4 — 0.1M© s _1 . In briefly, the abundances 
through the winds ejected from the O-rich and Si-rich disk layers are found 
to be not largely changed from those of the disk. On the other hand, for the 
winds from the He-rich and np-rich disk layers, the composition is largely 
altered from that of the disk. 

In Figure 4, we show the mass fractions of 56 Ni inside the winds from 
the accretion disk with M = 0.01M© s^ 1 . The abscissa is the radius from 
which the wind is launched. The abundances are evaluated at the time when 
T C j ~ 5 x 10 8 K. Abundant radioactives have not decayed significantly until 
the epoch. We find that the winds from the np-rich disk layer are abundant 
in 56 Ni, as suggested by the several authors. 3,8 We also find that 56 Ni is the 
most abundant in the ejecta from not only the np-rich layer but also the He- 
rich and Fe-rich layers and inner parts of the Si -rich layer. The smallness 
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Table 1 . Averaged abundances in the Fe and He layer of the disk before decay. 
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of 56 Ni in the ejecta from 10r g is attributed to efficient neutron capture 
on iron peak elements. The electron fraction, Y c , is initially 0.4830 in 
the ejecta. 56 Ni is consequently largely depleted through neutron capture, 
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Figure 4. Mass fraction of 56 Ni of winds ejected from various disk radii and various 
ejection velocities. 

while more neutron rich nuclei, such as 58 Ni, 60 Ni, and 64 Zn, are abundantly 
synthesized. On the other hand, in the ejecta launched from 20r g (Y e is 
0.4994 initially), neutron capture on iron peak elements is less efficient in 
the ejecta. The mass fraction of neutrons in the ejecta from 10r g is larger 
than that of protons by 0.034, which is too small to conduct r-process 
successfully. 



4.2. Chemical Composition of Winds 

We average abundances in winds ejected from the inner region (r < 50r s ) 
of the disk to estimate the yields through the winds from the disk. The 
averaging procedure is the same as in our previous study 13 but included 
with a weight of 1 — exp(— r O j/50r g ), (or 0, whichever is smaller than ), 
which means more massive ejection from smaller disk radius less than 50r 9 . 
Figure 5 shows the yields via the winds from the disk with M = 0.001, 
0.01, and 0.05M Q s _1 . The total mass of ejecta is set to be 1M . 3 The 
profiles of yields are similar for different M, but the ejected masses heavier 
than Cu are larger for higher M. This is due to the neutron richness of 
inner regions of the disk. It should be noted that appreciable amounts of 
63 Cu and 64 Zn can be produced through the winds from the disk without 
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Figure 5. Yields through the wind ejected from the disk with 0.001 (filled circles with 
solid line), 0.01 (circles with dashed line), and 0.05A/q s _1 (squares with dotted line) 
after decay when the total mass through the wind is IMq. 



neutron-rich regions (for O.OOIM0S 1 ). 



4.3. Chemical Composition of Winds from Neutron-rich 
Regions of Disks 

We investigate abundances in winds from neutron-rich regions of the disk 
with M > O.IMq s _1 . We use the same procedure as in §4.1 but with the 
larger nuclear reaction network which includes various neutron-rich nuclei. 
Figure 6 shows the abundances in the wind after decays for the represen- 
tative case: the wind launched from r e j = 10r g of the accretion disk with 
M = 0.1M Q s _1 . The electron fraction and the entropy per baryon are 
initially 0.30 and 13.1/cb, respectively The ejection velocities of the wind 
are taken to be O.OIwk (solid line) and O.OOIwk (dotted line). We find that 
neutron-rich nuclei are abundantly synthesized in the winds. More massive 
nuclei can be produced in the slower wind. This is because slower winds 
have enough time for nuclei to capture neutrons. Nuclei with mass num- 
bers ~ 120 and 180 are abundant for cases with w e j = O.OIwk and 0.001«k, 
respectively. For faster winds than O.Iwk, on the other hand, light elements 
are dominant. The abundances in the winds from neutron-rich regions of 
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Figure 6. Abundances in the wind ejected from r e j = 10r 9 of the disk with M = 
O.IMq s -1 . The ejection velocities of the wind are taken to be O.OI^k (solid line) and 
O.OOIuk (dotted line). 

the disk are sensitive to the ejection velocity as well as the electron fraction. 
Although the electron fractions are comparable to or less than that of the 
disk in our wind model, the fractions may increase to be larger than 0.5 
if the wind is driven via viscosity. 17 The abundances are likely to strongly 
depend on hydrodynamics of the wind. The yields from the neutron-rich 
disk therefore are still uncertain. 

5. Summary 

We have investigated nucleosynthesis inside the accretion disk associated 
with GRBs and inside winds launched from an inner region of the disk 
using the one-dimensional disk and wind models and the nuclear reaction 
network. The initial composition of accreting gas is taken to be that of an 
O-rich layer of a 20 Mq star before the core collapse. We have found that 
the disk consists of five layers characterized by dominant elements: 16 0, 
28 Si, 54 Fe (and 56 Ni), 4 He, and nucleons, and the individual layers shift 
inward with keeping the overall profiles of compositions as the accretion 
rate decreases. 56 Ni are abundantly ejected through the wind from the Fe- 
rich, He-rich and nucleon-rich disk layers with the electron fraction ~ 0.5. 
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In addition to iron group elements, heavier elements than Cu, in particular 
Cu and 64 Zn, arc massively produced via the wind. Various neutron- 
rich nuclei can be produced through the wind from neutron-rich regions of 
the disk in our simple wind model, though the estimated yields have large 
uncertainties. 
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